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We report on the results of a three-flavor osciUation analysis using Super-Kamiokande I atmo- 
spheric neutrino data, with the assumption of one mass scale dominance (Ami2=0). No significant 
flux change due to matter effect, which occurs when neutrinos propagate inside the Earth for ^137^0, 
has been seen either in a multi-GeV i^e-rich sample or in a !/^-rich sample. Both normal and 
inverted mass hierarchy hypotheses are tested and both are consistent with observation. Using 
Super-Kamiokande data only, 2-dimensional 90 % confidence allowed regions are obtained: mixing 
angles are constrained to sin^ ^13 < 0.14 and 0.37 < sin^ 623 < 0.65 for the normal mass hierarchy. 



Weaker constraints, sin 
hierarchy case. 



< 0.27 and 0.37 < sin 62s < 0.69, are obtained for the inverted mass 



PACS numbers: 14.60.Pq, 96.50.S- 



I. INTRODUCTION 

Recently a number of experiments have shown evidence 
for oscillations of atmospheric .1., J., ^, 4 , 5 , J.| , solar 0,13, 
reactor Q, and accelerator neutrinos p^. 

In the standard oscillation picture, the three neutrino 
flavor eigenstates are related to the mass eigenstates by 
a 3 X 3 unitary mixing matrix U: 



i=i 



Ua 



(1) 



In this picture, neutrino oscillations can be described by 
six parameters: two independent Arn^ij (Arnfj, ^77123), 
three mixing angles {9i2, ^23, ^13), and a CP-violating 
phase S. The mixing matrix U of Eq. ^ can be written 
as a product of three rotations, each described by one of 
the mixing angles: 



U = 




C12 S12 

-S12 C12 

1 



(2) 



where "s" represents sine of the mixing angle and "c" 
represents cosine. 

The "1-2" matrix describes solar mixing; the "2-3" ma- 
trix describes atmospheric neutrino mixing. The "1-3" 
mixing is known to be small; the best current limits on 
^13 come from the CHOOZ experiment nl|. 

As yet, the Super-Kamiokande atmospheric neutrino 
oscillation flts have been done within a two-flavor oscil- 
lation framework |5l Il2|. In this paper, we explore the 
atmospheric data in the context of a three-flavor analysis. 

Among the remaining problems in neutrino physics 
that can be answered by oscillation experiments are 
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whether ^13 is non-zero, and whether the hierarchy is nor- 
mal or inverted, i.e. whether Am^^ is positive or negative. 
At baselines and energies appropriate for atmospheric 
neutrinos, the signature of a non-zero ^13 is a matter- 
enhanced excess of upward-going electron- like events and 
possible additional small rate changes of upward-going 
muon-like events with respect to two-flavor v^ — )■ i^r tran- 
sition. The expected effects on electron-like (and muon- 
like) event rates differ for normal and inverted mass hier- 
archy cases because the matter effect and the cross sec- 
tion differ for !/„ and i/«. 



II. THREE NEUTRINO OSCILLATION WITH 
ONE MASS SCALE DOMINANT 

In general, neutrino oscillations are driven by differ- 



ences of squared masses, 



ml, 



ml,i 



We have adopted 



the so-called "one mass scale dominance" framework: 



« \m% 



(3) 



This approximation is supported by experimental ob- 
servations of solar, reactor, atmospheric and accelerator 
neutrino oscillations. The advantage of this framework is 
that the number of parameters involved in neutrino oscil- 
lations is reduced to three: two mixing angles (6'23,^i3) 
and one mass squared difference Am? , 



2 _ 



'1,2- 



(4) 



The ignored oscillation effects driven by the smaller mass 



difference Ar 



'12 



which might be observed 



for a neutrino energy of a few hundred MeV, are known 
to be greatly reduced in the case that the initial flavor 
flux ratio of v^j./ue is 2 and sin^ 26*23 is close to one [13j. 
For oscillations driven by the dominant Am?, a similar 
"screening" effect holds, but more weekly, because the 
flavor ratio starts deviating from 2 at 1 GeV and reaches 
~ 3 at 10 GeV. EquationOcan hold both for m.\ 2 << "13 
(normal mass hierarchy) and Wg << m^ 2 (inverted mass 
hierarchy). We present tests of both cases in this paper. 
In this framework, in the case of non-zero ^13 the 
neutrino transition and survival probabilities in vac- 
uum, valid for down-going atmospheric neutrinos, are 



expressed as 



1 - sin^ 26*13 sin^ 



1.27Am2L 
E 



.2. . 2„^ .2 /'1.27Am2L 
sm 023 sm 26/13 sm 



£; 



P(^m^'^m) = 1 



-4 cos ^lasin 023 (1 — cos 6*13 sin 6*23 
2 il.27Am^L 



E 



(5) 



where L is neutrino travel length in km from the neu- 
trino production point in the atmosphere, E is neutrino 
energy in GeV, and Am^ in eV^. In the limit of zero 
6*13, these equations reduce to pure v^ ^^ Vt two-flavor 
oscillation. Because Viv^ -^ Ve) hi Eq. is a func- 
tion of sin^ 26*13, both sin^ 6*13 ~ and sin^ ^13 ^ 1 can 
satisfy electron flavor disappearance constraints from re- 
actor neutrino experiments. However, sin ^13 ~ 1 is in- 
consistent with an observed large deficit of atmospheric 



because V{u^ 



i) in Eq. |S1 becomes ~ 1 for this 



case. Therefore, sin^ 6*13 > 0.5 is not discussed here. 

For neutrinos traversing the Earth, oscillation prob- 
ability is calculated taking into account Earth's mat- 
ter potential due to the forward scattering amp litude of 
charged current Ve and v^ interactions [Ij, [13, Hal ■ We 
adopted a model in which the Earth is well-approximated 
by four layers of a constant matter density (core 1: 
R<1221 km, p=13.0 g/cm^, core 2: 1221<R<3480 km, 
p=11.3 g/cm^, mantle: 3480<R<5701 km, p=5.0 g/cm^, 
surface: 5701<R<6371 km, p=3.3 g/cm^). The method 
of calculating matter oscillation probabilities in constant 
density is based on [iTj . 

In constant matter density, P(t'^ — > I'e) can be de- 
scribed by replacing the mixing angle ^13 and lS.m? in 
Eq. Elas 



P(l/^ -^ Ve) = P(j^e -^ V^) 

923 sin ZC13 



Sin2 023sin2 20f|sin2'^-^^^"''^^ 



E 



(6) 



and the effective mixing angle is 



sin-" 261 



M 
13 

Acc 



sin^ 26* 



13 



(cos20i3-Acc/Am2)2 
2V2Gi.iVeP 



sin^ 26*13 



(7) 



where Gf is the Fermi constant, A^e is the electron densi- 
ties in the medium and p is neutrino momentum J18i] . The 
matter potential term Acc bas the same absolute value, 
but opposite sign for neutrinos and anti-neutrinos. Due 
to the matter effect, the Mikheyev-Smirnov-Wolfenstein 
(MSW) resonance happens in that the v^^ -^ v^ os- 
cillation probability in Eq. and [7| becomes large at 




Ev (GeV) 

FIG. 1: Oscillation probability of v^ — > v^, (or Vy, — > v^ 
transition. For both figures, the horizontal axis shows neu- 
trino energy and the vertical axis shows the zenith angle of 
neutrino direction; cosGi/ = — 1 and cosSi^ = correspond 
to vertically upward and horizontal directions, respectively. 
Angles with cos Gi^ < —0.84 correspond to neutrinos passing 
through the earth core layers. Assumed oscillation parame- 
ters are (Am^ = 2.5 x 10"^ eV^ sin^ 6I23 = 0.5, sin^ 6*13 = 0.04 
). The top figure assumes neutrino oscillation in vacuum and 
the bottom figure takes into account Earth's matter effect. In 
the bottom figure, three high probability (> 40%) regions are 
shown which correspond to the MSW resonance at 3 GeV in 
the core layer, the MSW resonance at 7 GeV in the mantle 
layer, and the enhancement due to the core-mantle transition 
interference |2(l | at the energy between the two MSW regions. 



3 '^ 10 GeV neutrino energy even for the case of small 
^13 miMlIlllll- To demonstrate the behavior of v^ 
oscillations. Fig. ^ shows the transition probability of v^ 
to Vyi^ after traversing the Earth. Note that for normal 
mass hierarchy, the v^. flux is resonantly enhanced, and 
there is no enhancement for anti-neutrinos; the situation 
is reversed for the inverted mass hierarchy. In matter, the 
solar and the interference terms modify the v^, enhance- 
ment in the resonance region by less than 5%, justifying 
our assumption of Am^j = 0. 



III. DATA SAMPLE 

Super-Kamiokande is a 22.5 kt fiducial mass water 
Cherenkov detector located at a depth of 2700m water 
equivalent in the Kamioka mine, Gifu, Japan. The detec- 
tor is optically separated into two concentric cylindrical 
regions. The inner detector (ID) is instrumented with 
11,146 20-inch photomultipher tubes (PMT). The outer 
detector (OD) is instrumented with 1,885 8-inch PMTs. 
Details of the detector can be found in Ref . |23| . Phys- 
ical quantities associated with a neutrino event such as 
the interaction vertex, the number of Cherenkov rings, 
the direction of each ring, particle identification (PID), 
momentum, and number of Michel electrons are recon- 
structed by using hit timing and charge distributions of 
Cherenkov ring images recorded by PMTs on the ID wall. 

Atmospheric neutrino data are categorized into fully- 
contained (FC), partially-contained (PC), and upward- 
going muons (UPfi). In the FC events, all of their 
Cherenkov light is deposited in the ID. In the PC events, 
there is an exiting particle that deposits visible energy in 
the OD. Neutrino interactions in the rock below the de- 
tector produce UP^ events for which muons either stop in 
the detector, or pass through the detector. In the present 
analysis, we use 1489 live-days of FC, PC and 1646 days 
UP/x neutrino data taken from May 1996 through July 
2001 during the Super-Kamiokande I period. 

We divide the FC sample into sub-GeV and multi-GeV 
subsamples according to the visible energy as Eyis < 1.33 
GeV for sub-GeV, E^is > 1.33 GeV for multi-GeV where 
the visible energy of an event is the total energy assum- 
ing all Cherenkov light is from electromagnetic showers. 
The sample is also divided into single-ring and multi-ring 
by number of reconstructed Cherenkov rings, and into e- 
like and /x-like by PID of the most energetic ring. Sub- 
and multi-GeV, e- and fi-\ike events from the single-ring 
sample, and sub- and multi-GeV /i-like events from the 
multi-ring sample are used in the analysis, as for the 
i^fj. -^ Vt two-flavor oscillation analysis j3- We also di- 
vide the PC sample into 'OD stopping events' and 'OD 
through-going events' according to energy deposited in 
the OD 12]. Finally, the UP/x sample is divided into up- 
ward stopping muons (entering and stopping in the tank) 
and upward through-going muons (passing through the 
tank). 

As mentioned in the previous section, an excess of 
upward-going v^ and/or v^ in the several GeV region is 
expected for certain oscillation parameter sets. To im- 
prove the sensitivity to this case, a t/g-enriched sample is 
selected from the multi-GeV multi-ring e-like sample and 
is used in addition to the standard oscillation analysis 
samples. The z^e-enrichcd selection is based on a likeli- 
hood analysis using PID likelihood, momentum fraction 
of the most energetic ring, number of muon decay elec- 
trons, and distance between muon decay electron and 
primary neutrino interaction position. Figure |21 shows 
distributions of Monte Carlo (MC) events used to obtain 
the i^e-enriched likelihood function. The v^ + Ve charged 



current (CC) and i/^ CC -I- v^ CC -I- neutral current 
(NC) interactions are separately shown. These distribu- 
tions are divided into 5 energy regions; 1.33-2.5 GeV, 
2.5-5 GeV, 5-10 GeV, 10-20 GeV, and > 20 GeV. Each 
normalized signal (background) distributions forms a 
probability density function PDF^'j^{x{) {PDFl^j^{x{)), 
where i denotes the four observed variables and j is five 
energy bins. Then the i/e-enriched selection criterion 
for each event is defined as X)i=i{log(^^^i.j^(^j)) ^ 
\og{PDFP^{x^))} > 0. As is shown in Table H the 



TABLE I; Breakdown of multi-GeV multi-ring e-like events 
in the MC sample before and after a i^e-enrichment based on a 
likelihood analysis. For each interaction mode, the number of 
events normalized to 1489 live-days and the fraction are shown 
assuming pure u^ ^^ Vr two-flavor oscillation (Am'^ = 2.5 x 
10"^eV^ sin^ ^23 = 0.5, sin^ 6I13 = ). Survival efficiency for 
i/e 4- J^e CC events is 66.8%. 
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V^ + Ve CC 


v^ -\- Vfj, CC 


NC 


total 


e-like events 










no likelihood cut 


630.8 


242.2 


258.0 


1131.1 




55.8% 


21.4% 


22.8% 


100% 


likelihood cut 


421.3 


49.1 


86.8 


557.3 




75.6% 


8.8% 


15.6% 


100% 


efficiency 


66.8% 









J^e + i^e fraction in the multi-GeV multi-ring e-like sample 
is improved to 75.6% after the likelihood cut. A summary 
of the number of observed and expected FC multi-ring e- 
like events is shown in Table Hill in the appendix. 



IV. OSCILLATION ANALYSIS 

The oscillation analysis is performed by comparing 
data with MC equivalent to 100 years of detector ex- 
posure. The atmospheric neutrino flux calculation from 
|24l | and neutrino interaction model (NEUT) ^, .24I are 
used to simulate interactions with the nuclei of water, 
or in the case of upward muons, the nuclei of the rock 
surrounding the detector. We use a GEANT-based full 
detector simulation to generate the MC neutrino events. 

We employ a x^ test to perform three-flavor oscillation 
analysis. All events are divided into 37 momentum or en- 
ergy bins (lO-f-5 bins for FC single- and multi-ring e-like, 
8-1-4 bins for FC single- and multi-ring /i-like, 4-1-4 bins 
for OD stopping and OD through-going PC, and 1+1 bin 
for upward stopping and through-going muons) . The re- 
constructed quantities used for momentum binning for 
the various event classes are: the observed momentum of 
the charged lepton for FC single-ring e- and /i-like events 
(Piep), the sum of the energies of the observed rings (re- 
constructed particles) considering particle mass for FC 
multi-ring e-like events [Etot), and the sum of visible en- 
ergies of observed rings for FC multi-ring /i-like and PC 
events {Ems). Note that the binning of the energy scale 
differs from that of |a,[l3> ^-^^ events in the multi-GeV 
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FIG. 2: Livetime normalized MC distributions used in the i/e-enriched likelihood selection for multi-GeV multi-ring e-like 
events. Plots correspond to 1.33 < Evis < 5 GeV, 5 < Evia < 10 GeV, and Evia > 10 GeV from top to bottom. From left to 
right are shown PID likelihood for the most energetic ring (more negative means more electron-like), momentum fraction of 
the most energetic ring, number of Michel decay electrons, and square root of distance between decay electron and neutrino 
interaction vertex divided by visible energy of most energetic ring in (cm/MeV)^". Solid histograms are Ue + Ve CC, dashed 
are backgrounds (NC -|- v^ CC -I- i/^ CC). The Ve or u^ CC signals tend to give more electron- like PID and higher momentum 
fraction. In contrast, backgrounds tend to give more muon decay electrons and the longer decay electron distance due to 
energetic muons produced by i^^ or p^ CC interactions. 



energy range are divided more finely in order to obtain 
better sensitivity to Atti^. Each momentum bin is also 
divided into 10 bins equally spaced between cos 9 = — 1 
and COS0 = +1 (—1 < cos0 < for UP/x events), where 
cos 8 is the cosine of the zenith angle of the reconstructed 
particle direction. Table IIIII summarizes the number of 
observed and expected FC, PC and UP/x events for each 
bin. The total number of bins is 370. The number of 
events in each bin is compared with expectation and a 
X^ value is calculated according to a Poisson probability 
distribution defined by the following expression: 



X = 



370 

E 



^lK.pU+Y.fr^i]-Kbs\ 



+2NlAn 



N", 

obs 



Ns.A^ + ^=ifr^' 



43 / X 2 

e,; 



■E^ 



(8) 



iV'l 

obs 
exp 



Number of observed events in n-th bin 
Number of expected events in n-th bin 
i-th systematic error term 
Systematic error coefficient 
1 sigma value of systematic error 



where 



The expectation N^^^ is calculated using MC events 
corrected by oscillation probability. Here, systematic un- 
certainty factors explicitly multiply the N^ . We consid- 
ered 45 systematic error sources, which come from detec- 
tor calibration, neutrino flux, neutrino interactions and 
event selection. Most of them are in common with those 
listed in [5|, and additional systematic uncertainties re- 
lated to backgrounds of the e-like sample and upward- 
going muons, and sample normalization of e-like events 
are estimated as listed in TablellTl Among 45 errors, only 
43 contribute to the x^, because the overall normaliza- 
tion and sample normalization of multi-GeV multi-ring 
e-like are allowed to be free. The /" values are calculated 
and tabulated in advance for every bin and for every sys- 
tematic error source. A global scan is carried out on a 
(log]^o(^"^^)i sill ^23, sin 6*13) grid minimizing y^ with 
respect to 45 systematic error parameters. The e^ values 
are fit in order to minimize the x^ value. We used e^ 
such that the first derivative of x^ with respect to ei is 



TABLE II: Sources of systematic errors in addition to those 
in common with the z/^ <-> v^ oscillation analysis ,5]. Non-i/ 
backgrounds in UP^ samples are treated as fitting param- 
eters in this analysis, while they are taken into account by 
modifying statistical error size in Ref. Q- 

estimated error size (%) 
Non-i/ background " 

upward through-going muons 
upward stopping muons 
Non-(!/e CC) background ' 

multi-GeV single-ring e-like 
multi-GeV multi-ring e-like 
sample normalization 

of multi-GeV multi-ring e-like 
OD stopping PC/through-going PC 
separation 



3.0 

17 

14 
20 
free 

12 



"Cosmic ray muon backgrounds are assigned as systematic errors 
to the most horizontal zenith angle bins (—0.1 < cos < 0). 
''Contamination of Va CC interactions. 



grid point of (Am^, sin^6l23, sin^ 6113)= (-2.5xlO-3eV2, 
0.525, 0.00625). There is little difference in the x^min val- 
ues of the normal and inverted hierarchy cases; therefore 
both hypotheses are allowed by Super-K data. Figure 
Ul shows x^ — x^min distributions projected to sin 6*13, 
in which minimum x^ values for each Am^ and sin ^23 



are plotted. It is shown that the x 



X 



distribu- 



tion for the inverted hierarchy is flatter and a larger 
sin^ 013 value is allowed. The constraint on sin^ 6*13 is 
weaker for the inverted hierarchy case; sin^ 6*13 < 0.27 
and 0.37<sin^6'23<0.69 are allowed at 90% confidence 
level. 

The present analysis obtained upper limits on ^13 
which confirms CHOOZ experiment 11] (shown by Fig. 
|31 and EJ and Palo Verde [23. These limits are also 
consistent with the recent result by the K2K experi- 
ment [231 giving the upper hmit of sin^ 6'i3 ~ 0.06 at 
Am^ = 2.8 X 10"^eV^, assuming sin^ 6*23 = 0.5. 



zero (-#- = 0), which can be obtained by solving linear 
equations 26]. Since this equation has non-linear terms 
for our x^ definition, we use an approximate solution ob- 
tained by an iteration method. 

A global scan of the oscillation parameter grid as- 
suming normal mass hierarchy results in a minimum x^ 
value of x™™=377.39/368 DOF at the grid point (Am^, 
sin2 6l23, sin2 6'i3)= (2.5x10-3 eV^, 0.5, 0.0), which is 
consistent with f^ <-^- Vr two-flavor oscillation. Figure 
|31 shows the zenith angle distribution of each data sam- 
ple overlaid with non-oscillated and best-flt expectations. 
The fitted distributions agree well with data. Figure ^ 
shows the up-down asymmetry as a function of parti- 
cle momentum or total energy. The asymmetry (UP- 
DOWN)/(UP-hDOWN) distributions are consistent with 
the fitted expectation. No significant excess due to mat- 
ter effect is seen in the upward-going multi-GeV e-like 
sample, suggesting no evidence for non-zero ^13. Allowed 
regions of neutrino oscillation parameters are obtained 
based on the x^ defined in Eq. |H1 The 2-dimensional 
90 % (99 %) confidence level allowed region is defined 
to be X^=X^min+4.6 (9.2) and obtained as shown in 
Fig. El The region corresponding to sin^0i3<O.14 and 
0.37<sin^ 023<O.65 is allowed at 90% confidence level. 

Finally, we tested the inverted mass hierarchy hypoth- 
esis. Water Cherenkov detectors, such as Super-K, can- 
not discriminate between neutrinos and anti-neutrinos on 
an event-by-event basis. However, mass hierarchy affects 
the expected number of e-like events. Because of the 
lower cross section of Pe, an enhancement of the multi- 
GeV i/g-rich sample is expected to be suppressed and 
therefore the constraint on 6*13 will be weakened for the 
inverted mass hierarchy case. Allowed regions assum- 
ing inverted mass hierarchy are also obtained and shown 
in Fig. El x^„i„=377.31/368 DOF is obtained at the 



V. CONCLUSION 

In summary, a three-flavor oscillation analysis assum- 
ing one mass scale dominance (Ato^j =0) was per- 
formed with Super-Kamiokande I FC-l-PC-l-UP/z com- 
bined dataset. A multi-ring e-like sample, selected 
using a likelihood method, was newly introduced to 
increase the statistics of electron neutrinos and im- 
prove the sensitivity to ^13. The best-fitted parame- 
ters for three-flavor oscillation becomes (Am^, sin^6'23, 
sm^Bi^i)^ (2.5xl0~3 eV^, 0.5, 0.0) and the region of 
sin^ 6'i3<0.14 and 0.37<sin2 6'23<0.65 is allowed at 90 % 
confldence level, assuming normal mass hierarchy. We 
also tested the inverted mass hierarchy case: a wider re- 
gion, sin^ 6'i3<0.27 and 0.37<sin^ 6'23<0.69 is allowed at 
90 % confidence level. Both mass hierarchy hypotheses 
agree with our data. We obtained the upper limit on ^13, 
which is consistent with CHOOZ, Palo Verde, and K2K 
experiments, by using high statistics atmospheric neu- 
trinos. In contrast to these past experiments, the earth 
matter effect plays important role in the analysis. 
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A. Appendix 

Table IIIII summarizes the number of observed and ex- 
pected FC, PC and UP^ events for each bin. The Monte 
Carlo prediction does not include neutrino oscillations. 
These binned data are used in the oscillation analysis. 
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FIG. 3: Zenith angle distributions of FC e-like, /x-like, PC, 
and UP/i are shown for data (filled circles with statistical 
error bars), MC distributions without oscillation (boxes) and 
best-fit distributions (dashed). The non-oscillated MC shows 
the distribution without fitting and the box height shows the 
statistical error. In the case of non-zero &13, matter enhanced 
excess of electron-like events is expected in the zenith angle of 
— 1 < cos6' < —0.2 regions in the multi-GeV 1-ring and multi- 
ring electron-like samples. The i^^ in the resonance regions 
populate mainly in the multi-GeV single-ring muon, multi- 
ring muon, two PC, and UP stopping ^ samples. 
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FIG. 4: Asymmetry (UP-DOWN)/(UP+DOWN) as a func- 
tion of particle momentum for FC single-ring e-like events 
(top) and as a function of total energy for FC multi-ring e- 
like events (bottom), where UP (DOWN) refers to the num- 
ber of events in -1.0<cose<-0.2 (0.2<cose<1.0). Filled 
circles represent data (error bars are statistical), the line rep- 
resents best-fit distributions under a normal hierarchy as- 
sumption, and the filled area on the best-fit line represents 
the expected excess due to matter effect for (Am^, sin^6'23, 
sin^eia)^ (2.5x10-^ eV^ 0.5, 0.04). 
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FIG. 5; 90 % (thick line) and 99 % (thin hne) confidence level 
allowed regions are shown in sin'^ ^13 vs sin^ 6123 (top), Am^ 
vs sin^ 023 (middle), and Am^ vs sin^ S13 (bottom). Normal 
mass hierarchy {Am^>0) is assumed. The shaded area in 
the bottom figure shows the region excluded by the CHOOZ 
reactor neutrino experiment. 
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FIG. 6; 90 % (thick line) and 99 % (thin hne) confidence level 
allowed regions assuming inverted mass hierarchy (Ar7i^<0), 
shown in the same orientation as in Fig. |K| 
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Upward stopping muon 



28(51.2) I 23(54.1) | 37(56.7) | 30(65.0) | 27(67.6) | 37(68.2) | 37(78.9) | 48(81.0) | 65(94.0) | 85.7(96.9) | 

Upward through-going muon 
85(96.1) 1 113(115) 1 116(122) 1 138(137) 1 159(146) 1 183(169) 1 178(187) 1 267(211) 1 286(229) |316.6(257) | 



TABLE III: Summary of the number of observed (MC expected) FC, PC and UP/i events for each bin. Neutrino oscillation 
is not included in the MC prediction. Roman numbers I, II, ... X represent zenith angle regions —1 < cosO < —0.8, 
-0.8 < cose < -0.6, ... and 0.8 < cosG < 1.0 respectively for FC and PC events, -1 < cosG < -0.9, -0.9 < cosO < -0.8, 
... and 0.1 < cosB < 0.0 respectively for upward stopping and through-going muon events. The numbers 1 to 5 in the Piep 
column correspond to the momentum ranges <250, 250-400, 400-630, 630-1000 and >1000 MeV/c for sub-GeV samples and 
the numbers 6 to 10 correspond to <2.5, 2.5-5.0, 5.0-10, 10-20 and >20 GeV/c for multi-GeV samples. The letters a to f in 
the Etot and Evis columns correspond to energy ranges 0.2-1.33, 1.33-2.5, 2.5-5.0, 5.0-10, 10-20, >20 GeV. 



